Theory of superradiant scattering of laser light from Bose-Einstein condensates 
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In a recent MIT experiment, a new form of superradiant Rayleigh scattering was observed in Bose- 
Einstein condensates. We present a detailed theory of this phenomena in which condensate depeletion 
leads to mode competition, which, together with the directional dependence of the scattering rate, is 
ultimately responsible for the observed phenomena. The nonlinear response of the system is shown 
to be highly sensitive to initial quantum fluctuations which cause large run to run variations in the 
observed superradiant pulses. 
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With the recent advent of Bose-Einstein condensa- 
tion (BEC) in low-density alkali vapors j|§, a laser-like 
source of coherent monochromatic atomic matter-waves 
is now readily available. As the electromagnetic vacuum 
itself provides a nonlinear medium for atomic fields, an 
atomic BEC is thus an ideal system to study nonlinear 
wave-mixing and related phenomena. Indeed, nonlinear 
atom optics |j|-|5| is now an experimental reality with the 
recent observation of atomic four- wave mixing in conden- 
sate systems In addition to wave-mixing between 
atomic matter waves, the ability to generate laser-like 
atomic fields also raises the possibility to observe direct 
wave-mixing between atomic and optical fields. 

In the case where the atoms interact only with far 
off- resonant optical fields, the dominant atom-photon in- 
teraction is two-photon Rayleigh scattering. When the 
atoms are described as matter waves, Rayleigh scatter- 
ing is formally equivalent to a cubic nonlinearity, and 
therefore leads to four-wave mixing between atomic and 
optical fields. In recent work |g|-|lO[ along these lines, 
the scattering of light by a condensate from a strong 
pump laser into a weak quantized optical cavity mode 
was considered. This work was an extension of the collec- 
tive atomic recoil laser (CARL) 0,H into the regime of 
BEC, and focused on exploiting an instability in the light- 
matter interaction to parametrically amplify atomic and 
optical waves as well as to optically manipulate matter- 
wave coherence properties and generate entanglement be- 
tween atomic and optical fields. 

Recent experiments by Ketterle and coworkers ar MIT 
|u|, however, have demonstrated that this instability 
can play an important role also in the case in which 
laser light is scattered into the vacuum modes of the 
electromagnetic field. In these experiments, a variation 
of Dicke superradiance jl4| was observed in which the 
role of electronic coherence, which stores the memory of 
previous scattering events, is replaced by coherence be- 
tween center-of-mass momentum states, i.e. interference 
fringes in the atomic density. In this paper we present a 
multi-mode theory of condensate superradiance. Begin- 
ning with the elimination of the radiated light field as in 
the Wigner-Weiskopf theory of spontaneous emission, we 



then derive a linearized model which describes amplifi- 
cation of quantum fluctuations. This is then coupled to 
a 'classical' nonlinear model in which mode competition 
quenches scattering in all but the direction(s) of maxi- 
mum gain. The initial quantum fluctuations are shown 
to strongly influence superradiant pulse formation, and 
lead to large fluctuations between runs with identical ex- 
perimental parameters. In the MIT experiments this is 
clearly demonstrated by the presence of random spots in 
the angular distribution of the scattered photons. 

Our model consists of a Schrodinger field of two-level 
atoms coupled via the electric-dipole interaction to a far- 
off resonant pump laser field, as well as to the vacuum 
modes of the electromagnetic field. The pump laser has 
frequency ujq, wavevector ko = (wo/c)y, and its polariza- 
tion is taken along the x axis. Due to the large detun- 
ing between the pump frequency and the atomic tran- 
sition frequency u> a , we can eliminate the excited state 
field, and describe the atoms as a scalar field of ground 
state atoms. This atomic field is self-interacting, due to 
ground-ground collisions, however, we note that collisions 
which transfer populations are generally nonresonant and 
should make only a small contribution to the dynamics. 
The remaining collisions then simply give a mean field 
shift to the resonance frequency for quasiparticle exci- 
tations. As the effect of these shifts on our model are 
negligible, at present we include collisions only implicity 
in the determination of the condensate wavefunction. 

The effective Hamiltonian which describes the coupling 
of the atomic and electromagnetic fields is given by 

H = J d s r&(r)H (r)$(r) + J d 3 k^j(k)S t (k)S(k) 

d 3 kd 3 r [hg(k)W\r)tf(k)e i( - k °- k > r $(r) 
+ H.c] (1) 

where i'(r) is the atomic field operator, and 6(k) is the 
annihilation operator for a photon in mode k in the frame 
rotating at the pump frequency luq. The photon energy in 
this frame is given by k/(k) = c|k| — ujq. The single-atom 
Hamiltonian is given by H (r) = -(h 2 /2m)V 2 + V(r) + 
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/i|£!o| 2 /2A, V(r) being the trap potential, Qq the pump 
Rabi frequency, and A = ujq — uj a the pump detuning. 
The Hamiltonian (]]]) includes only scattering of pump 
photons, i.e. multiple scatterings between vacuum modes 
are neglected. The coupling coefficient for Rayleigh scat- 
tering between the pump and vacuum modes is 



.9(k) 
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where d is the magnitude of the atomic dipole moment 
for the transition involved. 

The atomic field is initially taken to be a number state 
in which N atoms occupy the trap ground state tpo(v), 
which satisfies [Ho(r) — Tifi] (fo(r) = 0, /i being the en- 
ergy of the trap ground state. The effect of atomic recoil 
during Rayleigh scattering between the pump and the 
vacuum mode k is therefore to transfer atoms into the 
state <fio(r) exp[i(ko — k) • r]. This suggest to expand the 
atomic field operator onto quasi-modes according to 



£(M)=$>|q)e- 



(3) 



where (r|q) = <po (r) exp (iq • r) , and w q = h\q\ 2 /2m. 
This is similar to the slowly varying envelope approx- 
imation from optical physics, the envelope being given 
here by <^o( r )- 

A discrete quantization of the q values follows from the 
requirement that the operators {c q } obey boson com- 
mutation relations [c q ,c q/ ] = (q|q'} ~ <5q,q'- Due to 
the finite size of the ground state wavefunction </?o(r), 
this means that q and q' must be seperated in k-space. 
Hence, the summation in Eq. (||) is taken to include the 
condensate mode q = as well as a grid of q-values as 
closely spaced as is consistent with orthogonally. Clearly, 
this expansion is not rigorously orthogonal and complete, 
however, it is sufficient to account for the quantum sta- 
tistical effects which occur above the critical phase-space 
density. 

An important aspect of BEC superradiance is the gen- 
eration of families of higher-order sidemodes due to the 
scattering of pump photons by the first-order sidemodes. 
For the scope of this paper, however, we consider a sim- 
plified model containing only the primary Rayleigh scat- 
tering process whereby a condensate atom is transferred 
to a first-order sidemode by scattering a pump photon. 
With this simplification, we insert the expansion (g) into 
Eq. (|l|) and arrive at the effective Hamiltonian 
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%(k)p q (k)e i ^ t c q &t(k)c + H.c] , (4) 



where p q (k) = J d 3 r |<^o(r)| 2 exp[— i(k — ko +q) • r] is the 
Fourier transform of the ground state density distribution 
centered at k = k — q. 



From the Hamiltonian (|4|) it is straightforward to de- 
rive the equation of motion for 6(k) , which upon formal 
integration yields 



b(M) = 6(k,0)e 



-iuj(k)t 



:£«7(k)p q (k)e* 



q#0 



x /*dTe- l( " (k)+ ^ )T c q (t-T)c (i-r), (5) 
Jo 

where the first term gives the free electromagnetic field, 
i.e. vacuum fluctuations, and the second term is the radi- 
ation field due to Rayleigh scattering. A nonzero expec- 
tation value of the coherence operator c q co indicates the 
presence of interference fringes, hence the radiated field 
increases as fringes build up. This term therefore leads 
to an instability where the memory of previous scattering 
events, stored in the matter-wave interference fringes, en- 
hances the present rate of Rayleigh scattering. 

Equation (^|) is then substituted into the equation of 
motion for c q .In the Markoff approximation, familiar 
from the Wigner-Weisskopf theory of spontaneous emis- 
sion, this yields 
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where 



G q = 2^y d 3 k|.g(k)| 2 |p q (k)| 2 <5Mk)+ W q) 



(6) 



(7) 



is the single-atom gain. In deriving Eq. (|6|) we have used 
the orthogonality of the states {|q}} to make the ap- 
proximation p*(k)p q '(k) ps | pq (k) | 2 <5q,q' , and neglected 
the principal part which accompanies the <5-function. If 
included, the principle part would contribute additional 
ground-state collisions due to the dipole-dipole interac- 
tion. 

For a closed atomic system, the total number of atoms 
is conserved, hence c £q — N — X^o^q^q- ^ or vei T 
short times we can therefore take cjco ~ N. In this case 
Eq. (0) reduces to 



dt 



= ^Nc a 



flit)- 



(8) 



where / q (t) is a noise operator whose correlation func- 
tions are given in the Markoff approximation by 



(4W/q(0> = 0, 
</qW/q(*')} - G^N6(t 



t>). 



(9) 



These noise operators allow the system to be triggered by 
quantum fluctuations, and hence describe "spontaneous" 
scattering which occurs in the absence of any sidemode 
population. 
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Equation (||) can be solved exactly, giving 
S q (t) = eW> N %(0) + f dre^/^fKt-r). (10) 



From Eq. (|10|) it is possible to compute the probability 
Pq(n,t) of having n atoms in mode q at time t, assum- 
ing zero population at t = 0. To accomplish this we first 
compute the anti-normally ordered characteristic func- 
tion, Xq.(v) = (exp[-?7*c q ]exp[?7c ( ! 1 ]), yielding 



Xq(r?) = e -lnl 2 (^( i )+D ! 



(11) 



where rj q (t) = exp(G q Ai) — 1 is the mean population 
of mode q at time t. We can identify expression Jll| ) 
as corresponding to a chaotic field The number 

distribution for a chaotic field is given by 



P q (n,i) 



,(*) 



»(*) 



-(n+l) 



(12) 



which for n q (t) > 1 is well approximated by 
exp[— n/n q (i)]/n q (f). 

When the mean population of a field mode is suffi- 
ciently large, correlation functions effectively factorize to 
all orders, and it becomes possible to formulate a 'clas- 
sical' description of the field dynamics. In the classi- 
cal theory, we can consider the sidemode populations 
as c-numbers and neglect the influence of the quantum 
noise operators as spontaneous scattering of atoms is now 
negligible in comparison to the stimulated contribution. 
With the inclusion of condensate depletion, the sidemode 
populations then obey the equations 



-n q = G n (N - 



q- 



(13) 



Momentum conservation tells us that for each atom scat- 
tered into the sidemode q, there is a photon scattered 
roughly in the direction k = ko — q. Hence, N x dlq/dt 
is the ideal photon count rate generated by the |0) — » |q) 
atomic center-of-mass transition. 

The 'classical' nonlinear model is applicable when 
n q (i) ^ 1, whereas the linearized quantum theory re- 
quires S q j ™qW ^ N. Provided that the number of 
active modes is not comparable to N, there is a signifi- 
cant overlap in the validity regimes of these two models. 
We then join them by choosing initial conditions for Eqs. 
( |l3| ) from Pq(n, t c ]), where t c j satisfies 1 -C « q (i c i) <SC N. 
Because the response is still linear at time t c j the result- 
ing nonlinear evolution does not depend on the particular 
choice of t c y 

We now analyze the geometrical dependence of the 
single-atom gain given by Eq. (0), and show that it is 
largest for radiation along the long-axis of the conden- 
sate. We first note that G q depends on g(k), which con- 
tains the dipole radiation pattern, as well as on p q (k), 
which depends on the geometry of the initial condensate. 



For a cigar-shaped condensate, aligned along the z-axis, 
p q (k) is a disc which lies parallel to the x-y-planc in k- 
space. The dimensions of the disc in k-space are roughly 
the inverse of the condensate dimensions in r-space. Thus 
for a condensate whose dimensions are large compared to 
an optical wavelength, the dimensions of p q (k) are small 
compared to ko. 

Since <?(k) is slowly varying compared to p q (k), it can 
be removed from the integral in Eq. (R), and evaluated 
at the center of p q (k). In addition we neglect the recoil 
shift w q in the ^-function as it has negligible effect on 
the value of G q . The remaining integral then defines the 
solid angle f2 q for the scattered radiation associated with 
the qth mode according to 
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d 3 k|p q (k)| 2 ( 5(|k|-A :o 
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which shows that only q values for which the center of 
p q (k) lies at a distance fco from the origin experience gain, 
a consequence of energy conservation. Thus for every 
active quasi-mode q there is a corresponding radiation 
direction k, such that q = fco(y — k). 

We can obtain a good estimate for fi q by taking 
|p q (k)| 2 to be an ellipsoid solid with the inverse dimen- 
sions of the condensate. This gives 



47T 

k 2 W 2 



cos 



k.z 



+ (L/W) 



2 sin 2 I 



k.z 



-1/2 



(15) 



where L is the length of the condensate along the z axis, 
W is the radial diameter, and 8^ - is the angle between 
the radiation direction and the long-axis of the conden- 
sate. Thus fi q is maximized for 0s. = 0,7r, correspond- 
ing to radiation along z and — z, where it is given by 
Sl q = 4:ir/k 2 W 2 . As 9^ - moves away from the z axis, f2 q 
is relatively flat until it reaches the geometric angle W/L, 
after which it drops off rapidly. It is important to note 
that for the isotropic case L = W there is no preferred 
direction, and a ring of radiation is instead observed. 

Taking into account all of these considerations, the ex- 
pression for the single-atom gain becomes 



Go 



sin 



cos z 



'k.z 



(L/wy 



(16) 



sin 



'k,z 



where Q = 3|rio| 2 r/8|A| 2 fc 2 IF 2 is the maximum single- 
atom gain, r = kqd 2 /3irheo being the single-atom spon- 
taneous decay rate, and 8^ ~ is the angle between the 
radiation and polarization directions. For the parame- 
ters of the MIT experiment we find Q ~ 4 x 10~ 4 • /, 
where / is the laser intensity in mW/cm 2 , and Q is given 
in Hz. A rough estimate of the duration of a superradi- 
ant pulse for the case = 10 6 and I = 100 mW/cm 2 
is t = \n(N)/QN ~ 150 /is, in excellent agreement with 
experimentally observed time scales. 

The interplay between the dependence of G q on the ra- 
diation direction and the nonlinearity in Eq. (O) leads 
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to mode competition, the outcome of which depends sen- 
sitively on the initial quantum fluctuations. When modes 
with different values of G q compete, the competition is 
'unfair' and the mode with the largest G q generally de- 
pletes all of the condensate atoms before the populations 
of the other modes have a chance to grow. Modes with 
the same G q , such as the quasi-modes corresponding to 
radiation along the z and — z directions, instead compete 
'fairly', and while the interplay is highly sensitive to the 
random initial conditions, neither mode necessarily wins, 
i.e. there is no real 'winner-takes-all' effect, and often a 
'tie' will occur. 

The angular dependence of the scattered light from a 
typical simulation is shown in Fig. 1, where we have plot- 
ted the photon count density as recorded by an ideal de- 
tector array located at a distance Z from the center of the 
condensate along the z axis. The dark regions correspond 
to maximum light intensity, while the white regions in- 
dicate negligible intensity. The center of the figure lies 
along the symmetry axis of the BEC, and half- width of 
the box, given by ZW/L, corresponds to the geometric 
radiation angle. The simulation was performed for a con- 
densate with TV = 10 6 atoms in a BEC with a width of 
10 fim and a length of 100 /im, in rough agreement with 
the MIT experiment. Only matter- wave quasimodes cor- 
responding to radiation within twice the geometric angle 
of +z and — z were included in the simulation. This set of 
576 modes, is more than sufficient, however, as a buildup 
of significant population in modes corresponding to radi- 
ation outside of the geometric angle was never observed. 
This can be attributed to the fact that these modes have 
a much smaller gain, and thus compete 'unfairly' with the 
'endfire' modes. Within 6^ ~ < W/L we see the results 
of the interplay between quantum fluctuations and 'fair' 
mode competition. The pattern of dark spots indicates 
multimode superradiance, and exhibits variation on the 
scale of Z/koW, corresponding to the solid angle of ra- 
diation for an endfire mode. The pattern, which arises 
from the amplification of quantum fluctuations, varies 
randomly from run to run, an effect which has been di- 
rectly observed experimentally. 

In conclusion, we note that the quasi-mode popula- 
tions will experience losses as the recoiling atoms even- 
tually propagate out of the condensate volume. The 
lifetime of the quasi-mode, however, is on the order of 
T q = mLq/h\q\, where L q is the length of the conden- 
sate along q. These losses tend to destroy the coherence 
between condensate and quasi-mode, which accounts for 
the observation of a threshold for superradiance in the 
MIT experiment: for insufficient laser power, the growth 
of matter- wave coherence cannot overcome the losses. As 
this threshold is very small, we have considered only the 
situation far above threshold, in which case the losses 
may be neglected. 

In general, the rate of matter-wave decoherence in su- 
perradiance or CARL type experiments is given by the 
ratio between the recoil velocity and the matter-wave 




FIG. 1. A typical numerical simulation of condensate su- 
perradiance. The scattered photon intensity is plotted as 
seen by a detector array located at a distance Z along the 
symmetry axis of the cigar-shaped BEC. The blacg regions 
correspond to maximum, and the white regions to negligible 
intensity, respectively. The width of the figure is 2ZW/L, 
corresponding to radiation within the geometric angle of the 
condensate. 

coherence length. As a BEC is maximally coherent, its 
coherence length is given, as above, by its spatial extent. 
For a noncondensed atomic cloud, however, it is instead 
given by the thermal DeBroglie wavelength, which is in 
general small compared to that of the BEC. As a result, 
the threshold for superradiance is significantly larger, for 
which reason superradiance was not observed above T c 
in the MIT experiment. We remark that in the case 
of the CARL, the presence of an optical cavity provides 
additional feedback, which can compensate for the lack 
of atomic coherence and allow for instability and gain. 

Lastly, we remark that the Hamiltonian (Q) describes 
the creation of correlated atom-photon pairs, and is 
therefore analogous to the optical parametric amplifier 
(OPA), which generates entangled two-photon states for 
a variety of applications, e.g. tests of Bells inequality, 
quantum cryptography, and quantum teleportation. It 
should therefore be possible to perform analogous exper- 
iments using entangled atom-photon states generated in 
a BEC superradiance experiment. 
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